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We present the path integral formulation of a broad class of generalized diffusion processes. Em- 
ploying the path integral we derive exact expressions for the path probability densities and joint 
probability distributions for the class of processes under consideration. We show that Continuous 
Time Random Walks (CTRWs) are included in our framework. A closed expression for the path 
probability distribution of CTRWs is found in terms of their waiting time distribution as the solution 
of a Dyson equation. As the formalism naturally includes the treatment of functionals a generalized 
Feynman-Kac formula is derived. 
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Introduction: In addition to Langevin equations and 
Fokker-Planck equations, the concept of path integrals 
plays a central role in the description of stochastic pro- 
cesses [l|, [|J . The path integral formulation is based on 
the specification of a probability measure assigned to 
each realization of the process and, hence, encodes the 
complete statistical information on the process. It was 
pioneered by Onsager and Machlup who derived an ex- 
pression for the probability density of a path of a linear 
Gaussian process [1] . Later their work has been extended 
to include nonlinear drift and diffusion coefficients as well 
as coloured noise 

Despite the success of models based on Brownian mo- 
tion and diffusion processes, over the last two decades it 
has become apparent that many dynamical systems in di- 
verse fields, ranging from biology to physics, cannot aptly 
be described within this framework [5|-|7| • Various gener- 
alizations of diffusion processes have been considered to 
account for such anomalous diffusion. Some of these gen- 
eralized models consist of an additional stochastic pro- 
cess. Examples include the Continuous Time Random 
Walk (CTRW) [8] which has been applied to model a va- 
riety of different systems @. The CTRW is defined as a 
proper random walk with an additional random process 
governing the waiting time between successive jumps. If 
the waiting time distribution is assumed to be scale-free 
the CTRW can be described by the celebrated fractional 
Fokker-Planck equation (For a review, we refer the reader 
to [6] ) . The corresponding Langevin equation can be for- 
mulated by application of the concept of subordination 
. A path integral formulation of CTRWs has not been 
established so far. Such a formulation, however, is of 
particular importance since CTRWs in general are non- 
Markovian processes and as such insufficiently described 
by single-point probability distributions 10]. Another 
example where stochastic dynamical processes are com- 
bined with an additional random process is the field of 



superstatistics fTTj ] . 

This Letter presents a path integral approach to diffu- 
sion processes which are additionally driven by a further 
stochastic process. This class of processes comprises the 
CTRW and thus fills the gap of the lacking path integral 
formulation for CTRW's. 

Path integral representation: The class of processes un- 
der consideration in this Letter is defined by the one- 
dimensional discrete Langevin equation 

qk+i=qk + rN(q k ) + VrDR k + akr k , (1) 

where the first three terms on the right hand side de- 
scribe a standard Langevin equation with time step r, 
drift N(q k ) and Gaussian random variables R k with am- 
plitude y/~D. The anomalous contribution stems from the 
last term where a k and r k both are random variables. 
The transition amplitude of this process can be obtained 
as 

p{qk+i\qk,a k ) = 

J dq^^iq k [ qk + 1 ^ qk ^ T N(q k )-OL k r k )-^§-q 2 k ^) 

where we assumed the stochastic force R k to be Gaus- 
sian distributed with zero-mean and vanishing correla- 
tion (R k Ri) = 5jk- The probability for a path starting 
at go is given by iteration of Eq. ([2J 

9k+i(qk+i]qk,a k ; ...,q ,a ) = (3) 

where Vq = Yl k=0 ^ and So(q,q~) is the Martin- Siggia- 
Rose action of a diffusion process with drift [T2| , written 
in discretized form, 

N tD , 

So(q,q) = i^2{q k (q k +i - qk - rN{q k )) +i—q 2 k } ■ (4) 

k=0 
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To further evaluate the expression for the path proba- 
bility, let us assume that the r^ are Gaussian random 
variables with zero mean and variance {r^ri) — rQSki- 
After averaging with respect to the r k we find 

gN+i(qN+v, qn, "jv; ••; go, ao) = 

i £ h=0 { « fc + 1 _ i k - T N (« fc ) ) + — ■ — ^r^- ql } 



Vqe 1 



EN 
k=0 



The probability density for a specific path 
fN+i(qN+V,-;qo) is then given by averaging with 
respect to the stochastic process au 

fN+i(qN+x; -;qo) = (6) 

Vqe s °^Z(i^f N ,..^q*) 
We have introduced the characteristic function 

= £\. ..;ao)e iE *=°'' fcQ *. (7) 

ctjv a a 

which is the characteristic function of the process ft = 
af . We have denoted the distribution of a = [a_/v, •■, cxq] 
by p(a). 

Propagators: The representation ^ for the path prob- 
ability serves as the starting point to determine propa- 
gators for the anomalous stochastic process defined by 
Eq.JT]). As usual, the propagator is obtained by inte- 
gration of all possible paths with appropriate boundary 
conditions 

g{q N ;q ,N) = / T>qf N (q N , .., q u q Q ) . (8) 



Let us consider a process consisting of two pure diffusions, 
i.e. we set N(qk) — 0. We obtain 



G(q;q ,N) 



p(a) 



(9) 



Observe that the quantity S(N) = J2k=o a 1 pl a Y s a k e Y 
role for the diffusion process. After introduction of the 

pdf Pn{S) = ^2 a p( a )S s j2 N ~ 1 a 2 we can rewr it° Eq.© 

as 



g(q N ,qo,N) = Y,PN(S) 



s=o 



g 2t(N D + SQ) 

v /2vrr(AL> + SQ) 



(10) 



In comparison to the regular diffusion process which takes 
place in the physical time Nt the second diffusion process 
can be interpreted to occur in the random time St whose 
distribution is determined by p(a). A corresponding re- 
sult for a continuous jump-diffusion process was obtained 
in UM with different methods. 



Multiple-time propagators can be obtained from the 
propagators G{qN, qo, N) in a straightforward manner 



00 00 



G(qN,qM,qo) — ^ ^ Pn-m,m(Sn — Sm>Sm) 

S n =0S m =Q 

x G{q N ;qM, N - M)g(q M ; q , M). (11) 



where Pn-m m(Sn — Sm',Sm) denotes the probability 
•(5) distribution of the variables Sn — Sm — J2k=M a k an d 



fe=M ' 

Sm = J2k=o a fc' which can be evaluated using the prob- 
ability distribution p(a). The generalization to n-time 
pdfs is straightforward and is given in terms of the pdf 
PN-k,k-l,...,m(SN - St; Sk - Si; ..S m ) and the product of 
the propagators Q(q N ;q M ,M — N). 

Examples for the a-process: All results obtained so far 
hold for a general process a. In the following we consider 
two specific examples for this process. 

The case where the aj. are i.i.d. random variables with 
a common pdf h(a) is in the realm of superstatisics [Tlj | . 

Renewal processes are obtained when the process a k is 
a binary string. If we denote the number of zeros between 
two successive ones by t and assume them to be i.i.d. 
random variables with a common pdf W(t), the process 
a is a renewal process. In this case the process described 
by Eq.JIJ is a CTRW with internal dynamics [U, 03] 
which for N(q k ) = D = includes the standard CTRW. 

Path integrals for CTRWs: The starting point of our 
treatment of CTRW path integrals are Eqs. (|6]) and (0 
and we need to characterize the process a. For the sake of 
simplicity, we restrict ourselves to processes with ao = 1, 
i.e. processes that start with an event. The aim is to 
represent the probability of a specific sequence p(a) in 
terms of the waiting time distribution, i.e. the number of 
zeros between two successive ones. Let Wi y j denote the 
probability to have i — j — 1 zeros between ones at j and 
i, i > j and let furthermore Wij = 0. It follows that the 
survival probability, i.e. the probability that a one at j 
is not followed by a further one till i is 1 — Wij. 

To proceed it is convenient to introduce the 
probability ffe(l, ct!fe_i, cx±, 1) of truncated strings, 
[1, 0, 0, .., 1, .., 0, 1] which end with the event at k. Ac- 
cording to the definition of renewal processes the trun- 
cated densities fulfill the relation 



k 

+ y^CT fc ;^(l, OLl-r, 1) , 



1=0 



where we have defined 



o~ki — 6 ah ^8 ak _ 1 fi..5 ai+1 fi8 au iWkd ■ 



(12) 



(13) 



The density of the truncated strings for k = N depends 
on the densities for k < N. Iterative application of (fl"2")) 
allows us to present the truncated density for k = N 
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according to 

u N (l,a k -i, ....,ai, 1) = (14) 

N N I 

cjvo + ^ &Nivm + o'Ni&ii'O'l'o + — 

where N > I > I' > ... > 0. This series can be summed, 
yielding 



JV 



(1, a k -i, ai,l) = - cr]^ l V w , 



(15) 



z=o 



where £ denotes the unit matrix and <r is a matrix with 
elements . Relation (fT5j) is valid for k < N and implies 
that the truncated densities are defined by a Dyson equa- 
tion whose explicit solution is given by matrix inversion 
of E-cr. 

The probability density p(a^, ao) of the renewal 
process then is determined from the truncated densities 
on the basis of the relationship 

N 

p(a N ,...,a ) = y^7^|(l,a;_ ll ..,a 1 ,l) , (16) 

1=0 

where we have defined 

INI = <W,l<5jV,i + (1 — W / Af;)^a N ,0^ajv-i,0--^Q i + i,O^Q(,l • 

(17) 

The explicit representation of the probability density is 
then just 



N N 



(a N , ...,a ) = ^2^2lNk[E - cr} ( kl &iq . 



(18) 



k=0 1=0 



Again, an expansion of the matrix [E — &] yields 
a representation of the probability density of a string 
[0, .., 1, ..1, ..1] in the number of one's contained in the 
string. Concluding, we state the characteristic function 
of the renewal process a which can be easily assessed on 
basis of Eq.fllSJl: 

N N 

Z(rj N ,..,r] Q ) = ^2^2lNk[E - a% x 'aio , (19) 



k=0 1=0 



where now 

Iki = (1- Wju)e*» 



(20) 



The combination of Eq. (0 with Eq. (|T9"f provides the de- 
sired path-integral representation for the class of CTRWs 
under consideration. 

Since this path-integral representation is a rather con- 
densed representation of the CTRW we proceed to de- 
rive a more transparent formulation. It is based on 
Eqs. (IT21 and (jT8"]l . It is clear, that the probability 



distribution /jv+i^at+i, (Ziv, -Qo) should have an anal- 
ogous representation. We introduce the abbreviations 
Ko(qN] — ; <Zi) for the path probability of the MSR action 
(U]) and V(qi+i; qi) for the short time propagator of the 
action with ai = 1 



V{qi+i]qi) 



y/2lTT(D + Q) 



(gi + i-gi-«(gj)r 

e 2r(D+Q) _ (21) 



The starting point is the definition of the truncated dis- 
tribution 

r]N(qN+i,qN,—,qo) = J~] fjy(l, ajy-i, ■., oti, 1) (22) 
x 9N+i(qN+i',qN, l; qn-u «jv-i; ■•; <?o, l) • 

Using the relation (fT2j) we arrive at the Dyson equation 
for the truncated distribution 



Vk(qk+i,Qk, -,qo) = Sk,oV{qi,q ) 

k 

+ ^2 v ilk+i,qk)W k iK (q k , ..q l+1 )rn(qi + i, ..q ) ■ 



(23) 



1=0 



Employing Eq. ((T6|) we can establish the connection be- 
tween the distributions rj and /: 



/jv+ifev+i, QN, ••) Qo) = 



(24) 



JV 



^(1 - W N i)K Q (q N+1 ,.., qi +1 )r)i(qi + i,..q ) . 



1=0 



The relations (|23[) and (|24[) yields a representation of the 
path-probability of the CTRW in terms of Kg and V. 
As we will show below, this representation can be the 
starting for determining expectation values of the path. 

It is possible to combine both equations, (|2"3"|) and (j2~4")l 
to obtain the representation of the joint probability dis- 
tribution as 

N 

fN+i(qN+i,qN, ••) 90 

1=0 

[Sik ~ V(qi + x,qi)Wi k K (q h ..q k+1 )}^ 1) V(qi,q )S kfi . 

Thereby, the Dyson equation (|2"3"|) was solved explicitly. 

Expansion of the inverse matrix then yields an expan- 
sion of the path probability in terms of the number of 
events. The first term /(°) yields the contribution of 
paths during which no event occurred after the first. 

fwliilN+i; -;qo) = (l-WN,o)K (q N+ i; ...;q 1 )V{qi',qa) , 

whereas the first order term contains the contribution 
from paths with one event 



N 



/jv+i (<?«+i ;---;<?o) = X^ 1 ~ w n,i)- 



(27) 



1=1 



K (q N+ i; qi+i)V(qi+i; qi)Wi fi Ko{qi; qi)V(qn q ) . 
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Higher order terms are obtained in the same way. Note 
the formal analogy to the self-energy corrections in quan- 
tum theory. In this case the corrections to the free propa- 
gator are given by the anomalous contributions from the 
a-process. By considering processes which start with an 
event, we imply non-equilibrium conditions. Relaxing 
this condition would lead to the case of aging CTRWs 

0, 

Generalized Fokker-Planck Equation and Generalized 
Feynman-Kac Formula: We want to conclude with the 
derivation of the generalized Fokker-Planck equation 
for the CTRW-case. Apparently, we can also derive 
Feynman-Kac formulas |16|. Let us consider the func- 
tional 



P N +i(p, qw+i) 



(28) 



r N 

/ H dq^^o v(^f N+1 ( qN+1 , .., qo ) 

J 3=0 



The quantity Pjv+i(0, qN+i) is just the probability dis- 
tribution fN + i(qN + i) which obeys a generalized Fokker- 
Planck equation. Defining now 



G(q N+ i,qi;p,N,l) 



(29) 



r N 

J LI d qj e l P^ U ^f N+1 ^(q N+1 ,.., qi ) 



3=1+1 

as well as 



c»(p,?i+i)= n^^sfe.-)*). (3°) 

J 3=0 

we obtain from from Eqs. (|2U)) and (j2"4")l the relations 

N 

PN+i(p,qN+i) = J2( 1 - W »^ ( 31 ) 

k=0 

x / dq k +iG(q N+ i,q k+ i;p,N,k)Ci(p,qk+i) , 



k 

Cfc(P) gfc+l) = V(qi, 9o)Sl,o + ^ / dq k I dqi +1 (32) 
V(q k+1 ,q k )W ki iG(q k ,qi+i,k,l + 1)0 (p, qi+i) ■ 



Now Eqs. pip and (j52")l are just the time-discrete ver- 
sions of the equations which served as the starting point 
for the derivation of generalized Fokker-Planck equation 
for the CTRW with internal dynamics [l3[ eventually 
leading to 



P(p,q,t)= f dt'Qit-t'^e^'^P^q^'), 
Jo 

(33) 

where % = C + ipU(q). C is the generator of the 
Fokker-Planck process with a = and Q(t — t') is the 
common time-evolution kernel of CTRWs [6]. For the 
case of CTRWs without internal dynamics, i.e. £ = 0, 
this equation was recently derived in [T^ |. 

Conclusions and Outlook In summary, we have pre- 
sented the path integral formulation for a class of gen- 
eralized diffusion processes which includes the CTRW. 
Applying the path integral, we have derived general ex- 
pressions for the propagators and the multipoint distri- 
butions. For the case of CTRW-processes we found a 
closed form expression for the probability density of a 
path in terms of the waiting time distribution. We ex- 
pect that this expression will help to evaluate interesting 
functionals for anomalous diffusion processes. Further- 
more, we want to mention that it is tempting to consider 
the case, where the renewal process a is generated by a 
one-dimensional map [l8j . along similar lines. 
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